Proteomic Analyses in Diverse Populations Improved Risk Prediction and Identified New Drug Targets for Type 2 Diabetes.

OBJECTIVE: Integrated analyses of plasma proteomics and genetic data in prospective studies can help assess the causal relevance of proteins, improve risk prediction, and discover novel protein drug targets for type 2 diabetes (T2D). RESEARCH DESIGN AND METHODS: We measured plasma levels of 2,923 proteins using Olink Explore among ∼2,000 randomly selected participants from China Kadoorie Biobank (CKB) without prior diabetes at baseline. Cox regression assessed associations of individual protein with incident T2D (n = 92 cases). Proteomic-based risk models were developed with discrimination, calibration, reclassification assessed using area under the curve (AUC), calibration plots, and net reclassification index (NRI), respectively. Two-sample Mendelian randomization (MR) analyses using cis-protein quantitative trait loci identified in a genome-wide association study of CKB and UK Biobank for specific proteins were conducted to assess their causal relevance for T2D, along with colocalization analyses to examine shared causal variants between proteins and T2D. RESULTS: Overall, 33 proteins were significantly associated (false discovery rate < 0.05) with risk of incident T2D, including IGFBP1, GHR, and amylase. The addition of these 33 proteins to a conventional risk prediction model improved AUC from 0.77 (0.73-0.82) to 0.88 (0.85-0.91) and NRI by 38%, with predicted risks well calibrated with observed risks. MR analyses provided support for the causal relevance for T2D of ENTR1, LPL, and PON3, with replication of ENTR1 and LPL in Europeans using different genetic instruments. Moreover, colocalization analyses showed strong evidence (pH4 > 0.6) of shared genetic variants of LPL and PON3 with T2D. CONCLUSIONS: Proteomic analyses in Chinese adults identified novel associations of multiple proteins with T2D with strong genetic evidence supporting their causal relevance and potential as novel drug targets for prevention and treatment of T2D.

Plasma Proteomics to Identify Drug Targets for Ischemic Heart Disease.

BACKGROUND: Integrated analyses of plasma proteomic and genetic markers in prospective studies can clarify the causal relevance of proteins and discover novel targets for ischemic heart disease (IHD) and other diseases. OBJECTIVES: The purpose of this study was to examine associations of proteomics and genetics data with IHD in population studies to discover novel preventive treatments. METHODS: We conducted a nested case-cohort study in the China Kadoorie Biobank (CKB) involving 1,971 incident IHD cases and 2,001 subcohort participants who were genotyped and free of prior cardiovascular disease. We measured 1,463 proteins in the stored baseline samples using the OLINK EXPLORE panel. Cox regression yielded adjusted HRs for IHD associated with individual proteins after accounting for multiple testing. Moreover, cis-protein quantitative loci (pQTLs) identified for proteins in genome-wide association studies of CKB and of UK Biobank were used as instrumental variables in separate 2-sample Mendelian randomization (MR) studies involving global CARDIOGRAM+C4D consortium (210,842 IHD cases and 1,378,170 controls). RESULTS: Overall 361 proteins were significantly associated at false discovery rate <0.05 with risk of IHD (349 positively, 12 inversely) in CKB, including N-terminal prohormone of brain natriuretic peptide and proprotein convertase subtilisin/kexin type 9. Of these 361 proteins, 212 had cis-pQTLs in CKB, and MR analyses of 198 variants in CARDIOGRAM+C4D identified 13 proteins that showed potentially causal associations with IHD. Independent MR analyses of 307 cis-pQTLs identified in Europeans replicated associations for 4 proteins (FURIN, proteinase-activated receptor-1, Asialoglycoprotein receptor-1, and matrix metalloproteinase-3). Further downstream analyses showed that FURIN, which is highly expressed in endothelial cells, is a potential novel target and matrix metalloproteinase-3 a potential repurposing target for IHD. CONCLUSIONS: Integrated analyses of proteomic and genetic data in Chinese and European adults provided causal support for FURIN and multiple other proteins as potential novel drug targets for treatment of IHD.

Conventional and genetic associations of adiposity with 1463 proteins in relatively lean Chinese adults.

Adiposity is associated with multiple diseases and traits, but little is known about the causal relevance and mechanisms underlying these associations. Large-scale proteomic profiling, especially when integrated with genetic data, can clarify mechanisms linking adiposity with disease outcomes. We examined the associations of adiposity with plasma levels of 1463 proteins in 3977 Chinese adults, using measured and genetically-instrumented BMI. We further used two-sample bi-directional MR analyses to assess if certain proteins influenced adiposity, along with other (e.g. enrichment) analyses to clarify possible mechanisms underlying the observed associations. Overall, the mean (SD) baseline BMI was 23.9 (3.3) kg/m2, with only 6% being obese (i.e. BMI ≥ 30 kg/m2). Measured and genetically-instrumented BMI was significantly associated at FDR < 0.05 with levels of 1096 (positive/inverse: 826/270) and 307 (positive/inverse: 270/37) proteins, respectively, with FABP4, LEP, IL1RN, LSP1, GOLM2, TNFRSF6B, and ADAMTS15 showing the strongest positive and PON3, NCAN, LEPR, IGFBP2 and MOG showing the strongest inverse genetic associations. These associations were largely linear, in adiposity-to-protein direction, and replicated (> 90%) in Europeans of UKB (mean BMI 27.4 kg/m2). Enrichment analyses of the top > 50 BMI-associated proteins demonstrated their involvement in atherosclerosis, lipid metabolism, tumour progression and inflammation. Two-sample bi-directional MR analyses using cis-pQTLs identified in CKB GWAS found eight proteins (ITIH3, LRP11, SCAMP3, NUDT5, OGN, EFEMP1, TXNDC15, PRDX6) significantly affect levels of BMI, with NUDT5 also showing bi-directional association. The findings among relatively lean Chinese adults identified novel pathways by which adiposity may increase disease risks and novel potential targets for treatment of obesity and obesity-related diseases.

Multi-ancestry genome-wide association analyses improve resolution of genes and pathways influencing lung function and chronic obstructive pulmonary disease risk.

Lung-function impairment underlies chronic obstructive pulmonary disease (COPD) and predicts mortality. In the largest multi-ancestry genome-wide association meta-analysis of lung function to date, comprising 580,869 participants, we identified 1,020 independent association signals implicating 559 genes supported by ≥2 criteria from a systematic variant-to-gene mapping framework. These genes were enriched in 29 pathways. Individual variants showed heterogeneity across ancestries, age and smoking groups, and collectively as a genetic risk score showed strong association with COPD across ancestry groups. We undertook phenome-wide association studies for selected associated variants as well as trait and pathway-specific genetic risk scores to infer possible consequences of intervening in pathways underlying lung function. We highlight new putative causal variants, genes, proteins and pathways, including those targeted by existing drugs. These findings bring us closer to understanding the mechanisms underlying lung function and COPD, and should inform functional genomics experiments and potentially future COPD therapies.

Education interacts with genetic variants near GJD2, RBFOX1, LAMA2, KCNQ5 and LRRC4C to confer susceptibility to myopia.

Myopia most often develops during school age, with the highest incidence in countries with intensive education systems. Interactions between genetic variants and educational exposure are hypothesized to confer susceptibility to myopia, but few such interactions have been identified. Here, we aimed to identify genetic variants that interact with education level to confer susceptibility to myopia. Two groups of unrelated participants of European ancestry from UK Biobank were studied. A 'Stage-I' sample of 88,334 participants whose refractive error (avMSE) was measured by autorefraction and a 'Stage-II' sample of 252,838 participants who self-reported their age-of-onset of spectacle wear (AOSW) but who did not undergo autorefraction. Genetic variants were prioritized via a 2-step screening process in the Stage-I sample: Step 1 was a genome-wide association study for avMSE; Step 2 was a variance heterogeneity analysis for avMSE. Genotype-by-education interaction tests were performed in the Stage-II sample, with University education coded as a binary exposure. On average, participants were 58 years-old and left full-time education when they were 18 years-old; 35% reported University level education. The 2-step screening strategy in the Stage-I sample prioritized 25 genetic variants (GWAS P < 1e-04; variance heterogeneity P < 5e-05). In the Stage-II sample, 19 of the 25 (76%) genetic variants demonstrated evidence of variance heterogeneity, suggesting the majority were true positives. Five genetic variants located near GJD2, RBFOX1, LAMA2, KCNQ5 and LRRC4C had evidence of a genotype-by-education interaction in the Stage-II sample (P < 0.002) and consistent evidence of a genotype-by-education interaction in the Stage-I sample. For all 5 variants, University-level education was associated with an increased effect of the risk allele. In this cohort, additional years of education were associated with an enhanced effect of genetic variants that have roles including axon guidance and the development of neuronal synapses and neural circuits.

Non-additive (dominance) effects of genetic variants associated with refractive error and myopia.

Genome-wide association studies (GWAS) have revealed that the genetic contribution to certain complex diseases is well-described by Fisher's infinitesimal model in which a vast number of polymorphisms each confer a small effect. Under Fisher's model, variants have additive effects both across loci and within loci. However, the latter assumption is at odds with the common observation of dominant or recessive rare alleles responsible for monogenic disorders. Here, we searched for evidence of non-additive (dominant or recessive) effects for GWAS variants known to confer susceptibility to the highly heritable quantitative trait, refractive error. Of 146 GWAS variants examined in a discovery sample of 228,423 individuals whose refractive error phenotype was inferred from their age-of-onset of spectacle wear, only 8 had even nominal evidence (p < 0.05) of non-additive effects. In a replication sample of 73,577 individuals who underwent direct assessment of refractive error, 1 of these 8 variants had robust independent evidence of non-additive effects (rs7829127 within ZMAT4, p = 4.76E-05) while a further 2 had suggestive evidence (rs35337422 in RD3L, p = 7.21E-03 and rs12193446 in LAMA2, p = 2.57E-02). Accounting for non-additive effects had minimal impact on the accuracy of a polygenic risk score for refractive error (R2 = 6.04% vs. 6.01%). Our findings demonstrate that very few GWAS variants for refractive error show evidence of a departure from an additive mode of action and that accounting for non-additive risk variants offers little scope to improve the accuracy of polygenic risk scores for myopia.

Genome-wide association meta-analysis of corneal curvature identifies novel loci and shared genetic influences across axial length and refractive error.

Corneal curvature, a highly heritable trait, is a key clinical endophenotype for myopia - a major cause of visual impairment and blindness in the world. Here we present a trans-ethnic meta-analysis of corneal curvature GWAS in 44,042 individuals of Caucasian and Asian with replication in 88,218 UK Biobank data. We identified 47 loci (of which 26 are novel), with population-specific signals as well as shared signals across ethnicities. Some identified variants showed precise scaling in corneal curvature and eye elongation (i.e. axial length) to maintain eyes in emmetropia (i.e. HDAC11/FBLN2 rs2630445, RBP3 rs11204213); others exhibited association with myopia with little pleiotropic effects on eye elongation. Implicated genes are involved in extracellular matrix organization, developmental process for body and eye, connective tissue cartilage and glycosylation protein activities. Our study provides insights into population-specific novel genes for corneal curvature, and their pleiotropic effect in regulating eye size or conferring susceptibility to myopia.

Evidence That Emmetropization Buffers Against Both Genetic and Environmental Risk Factors for Myopia.

Purpose: To test the hypothesis that emmetropization buffers against genetic and environmental risk factors for myopia by investigating whether risk factor effect sizes vary depending on children's position in the refractive error distribution. Methods: Refractive error was assessed in participants from two birth cohorts: Avon Longitudinal Study of Parents and Children (ALSPAC) (noncycloplegic autorefraction) and Generation R (cycloplegic autorefraction). A genetic risk score for myopia was calculated from genotypes at 146 loci. Time spent reading, time outdoors, and parental myopia were ascertained from parent-completed questionnaires. Risk factors were coded as binary variables (0 = low, 1 = high risk). Associations between refractive error and each risk factor were estimated using either ordinary least squares (OLS) regression or quantile regression. Results: Quantile regression: effects associated with all risk factors (genetic risk, parental myopia, high time spent reading, low time outdoors) were larger for children in the extremes of the refractive error distribution than for emmetropes and low ametropes in the center of the distribution. For example, the effect associated with having a myopic parent for children in quantile 0.05 vs. 0.50 was as follows: ALSPAC: age 15, -1.19 D (95% CI -1.75 to -0.63) vs. -0.13 D (-0.19 to -0.06), P = 0.001; Generation R: age 9, -1.31 D (-1.80 to -0.82) vs. -0.19 D (-0.26 to -0.11), P < 0.001. Effect sizes for OLS regression were intermediate to those for quantiles 0.05 and 0.50. Conclusions: Risk factors for myopia were associated with much larger effects in children in the extremes of the refractive error distribution, providing indirect evidence that emmetropization buffers against both genetic and environmental risk factors.

Quantile regression analysis reveals widespread evidence for gene-environment or gene-gene interactions in myopia development.

A genetic contribution to refractive error has been confirmed by the discovery of more than 150 associated variants in genome-wide association studies (GWAS). Environmental factors such as education and time outdoors also demonstrate strong associations. Currently however, the extent of gene-environment or gene-gene interactions in myopia is unknown. We tested the hypothesis that refractive error-associated variants exhibit effect size heterogeneity, a hallmark feature of genetic interactions. Of 146 variants tested, evidence of non-uniform, non-linear effects were observed for 66 (45%) at Bonferroni-corrected significance (P < 1.1 × 10-4) and 128 (88%) at nominal significance (P < 0.05). LAMA2 variant rs12193446, for example, had an effect size varying from -0.20 diopters (95% CI -0.18 to -0.23) to -0.89 diopters (95% CI -0.71 to -1.07) in different individuals. SNP effects were strongest at the phenotype extremes and weaker in emmetropes. A parsimonious explanation for these findings is that gene-environment or gene-gene interactions in myopia are pervasive.